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Abstract. We show how fundamental ideas from signal processing, multiscale 
theory and wavelets may be applied to nonlinear dynamics. 

The problems from dynamics include iterated function systems (IFS), 
dynamical systems based on substitution such as the discrete systems built 
on rational functions of one complex variable and the corresponding Julia 
sets, and state spaces of subshifts in symbolic dynamics. Our paper serves to 
motivate and survey our recent results in this general area. Hence we leave 
out some proofs, but instead add a number of intuitive ideas which we hope 
will make the subject more accessible to researchers in operator theory and 
systems theory. 

Mathematics Subject Classification (2000). 42C40, 42A16, 43A65, 42A65. 
Keywords. Nonlinear dynamics, martingale, Hilbert space, wavelet. 



1. Introduction 

In the past twenty years there has been a great amount of interest in the theory 
of wavelets, motivated by applications to various fields such as signal processing, 
data compression, tomography, and subdivision algorithms for graphics (Our latest 
check on the word "wavelet" in Google turned up over one million and a half 
results, 1,590,000 to be exact). It is enough here to mention two outstanding 
successes of the theory: JPEG 2000, the new standard in image compression, and 
the WSQ (wavelet scalar quantization) method which is used now by the FBI 
to store its fingerprint database. As a mathematical subject, wavelet theory has 
found points of interaction with functional and harmonic analysis, operator theory, 
ergodic theory and probability, numerical analysis and differential equations. With 
the explosion of information due to the expansion of the internet, there came a 
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need for faster algorithms and better compression rates. These have motivated the 
research for new examples of wavelets and new wavelet theories. 

Recent developments in wavelet analysis have brought together ideas from 
engineering, from computational mathematics, as well as fundamentals from rep- 
resentation theory. This paper has two aims: One to stress the interconnections, 
as opposed to one aspect of this in isolation; and secondly to show that the funda- 
mental Hilbert space ideas from the linear theory in fact adapt to a quite wide class 
of nonlinear problems. This class includes random-walk models based on martin- 
gales. As we show, the theory is flexible enough to allow the adaptation of pyramid 
algorithms to computations in the dynamics of substitution theory as it is used 
in discrete dynamics; e.g., the complex dynamics of Julia sets, and of substitution 
dynamics. 

Our subject draws on ideas from a variety of directions. Of these direc- 
tions, we single out quadrature-mirror filters from signal/image processing. High- 
pass/low-pass signal processing algorithms have now been adopted by pure mathe- 
maticians, although they historically first were intended for speech signals, see 
Perhaps unexpectedly, essentially the same quadrature relations were rediscovered 
in operator-algebra theory, and they are now used in relatively painless construc- 
tions of varieties of wavelet bases. The connection to signal processing is rarely 
stressed in the math literature. Yet, the flow of ideas between signal processing 
and wavelet math is a success story that deserves to be told. Thus, mathematicians 
have borrowed from engineers; and the engineers may be happy to know that what 
they do is used in mathematics. 

Our presentation serves simultaneously to motivate and to survey a number 
of recent results in this general area. Hence we leave out some proofs, but instead 
we add a number of intuitive ideas which we hope will make the subject more 
accessible to researchers in operator theory and systems theory. Our theorems 
with full proofs will appear elsewhere. An independent aim of our paper is to 
point out several new and open problems in nonlinear dynamics which could likely 
be attacked with the general multiscale methods that are the focus of our paper. 

In Section [2] below, we present background material from signal processing 
and from wavelets in a form which we hope will be accessible to operator theorists, 
and to researchers in systems theory. This is followed in Section|3by a presentation 
of some motivating examples from nonlinear dynamics. They are presented such as 
to allow the construction of appropriate Hilbert spaces which encode the multiscale 
structure. Starting with a state space X from dynamics, our main tool in the 
construction of a multiscale Hilbert space H(X) is the theory of random walk and 
martingales. The main part of our paper, Sections 01 and [S] serves to present our 
recent new and joint results. 
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2. Connection to signal processing and wavelets 

We will use the term "filter" in the sense of signal processing. In the simplest 
case, time is discrete, and a signal is just a numerical sequence. We may be acting 
on a space of signals (sequences) using the operation of Cauchy product; and the 
operations of down-sampling and up-sampling. This viewpoint is standard in the 
engineering literature, and is reviewed in ^U] (see also JU| an d H2|) for the benefit 
of mathematicians. 

A numerical sequence (ajt) represents a filter, but it is convenient, at the 
same time, also to work with the frequency-response function. By this we mean 
simply the Fourier series (see (I3.3[) below) corresponding to the sequence (afc). 
This Fourier series is called the filter function, and in one variable we view it 
as a function on the one-torus. (This is using the usual identification of periodic 
functions on the line with functions on the torus.) The advantage of this dual 
approach is that Cauchy product of sequences then becomes pointwise product of 
functions. 

We will have occasion to also work with several dimensions d, and then our 
filter function represents a function on the d-torus. While signal processing algo- 
rithms are old, their use in wavelets is of more recent vintage. From the theory of 
wavelets, we learn that filters are used in the first step of a wavelet construction, 
the so-called multiresolution approach. In this approach, the problem is to con- 
struct a function on K or on K rf which satisfies a certain renormalization identity, 
called the scaling equation, and we recall this equation in two versions (|2.2|) and 
(|3.2I) below. The numbers (a^) entering into the scaling equation turn out to be the 
very same numbers the signal processing engineers discovered as quadrature-mirror 
filters. 

A class of wavelet bases are determined by filter functions. We establish a 
duality structure for the category of filter functions and classes of representations 
of a certain C*-algebra. In the process, we find new representations, and we prove 
a correspondence which serves as a computational device. Computations are done 
in the sequence space £ 2 , while wavelet functions are in the Hilbert space £ 2 (IR), 
or some other space of functions on a continuum. The trick is that we choose 
a subspace of L 2 (R) in which computations are done. The choice of subspace 
is dictated by practical issues such as resolution in an image, or refinement of 
frequency bands. 

We consider non-abelian algebras containing canonical maximal abelian sub- 
algebras, i.e., C(X) where X is the Gelfand space, and the representations define 
measures jionl. Moreover, in the examples we study, it turns out that X is an 
affine iterated function system (IFS), of course depending on the representation. 
In the standard wavelet case, X may be taken to be the unit interval. Following 
Wickerhauser et al. [HOj, these measures fi are used in the characterization and 
analysis of wavelet packets. 

Orthogonal wavelets, or wavelet frames, for L 2 (W. d ) are associated with 
quadrature- mirror filters (QMF), a set of complex numbers which relate the dyadic 
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scaling of functions on M. d to the Z d -translates. In the paper [531, we show that 
generically, the data in the QMF-systems of wavelets are minimal, in the sense 
that it cannot be nontrivially reduced. The minimality property is given a geo- 
metric formulation in I 2 ; and it is then shown that minimality corresponds 
to irreducibility of a wavelet representation of the Cuntz algebra On- Our result is 
that this family of representations of Oat on £ 2 (Z d ) is irreducible for a generic set 
of values of the parameters which label the wavelet representations. Since MRA- 
wavelets correspond to representations of the Cuntz algebras On, we then get, as 
a bonus, results about these representations. 

Definition 2.1. (Ruelle's operators.) Let (X, fi) be a finite Borel measure space, and 
let r : X — > X be a finite-to-one mapping of X onto itself. (The measure fi will be 
introduced in Section^bclow, and it will be designed so as to have a certain strong 
invariance property which generalizes the familiar property of Haar measure in the 
context of compact groups.) Let V: X — > [0, oo) be a given measurable function. 
We then define an associated operator R = Ry, depending on both V and the 
endomorphism r, by the following formula 

Rvf(x) = #r - 1(!c) E v(y)f(y), f&L\x,»). (2.1) 

r(y)=x 

Each of the operators Ry will be referred to as a Ruelle operator, or a transition 
operator; and each R v clearly maps positive functions to positive functions. (When 
we say "positive" we do not mean "strictly positive", but merely "non- negative".) 

We refer to our published papers/monographs |54l 1401 |4"5] about the spectral 
picture of a positive transition operator R (also called the Ruelle operator, or the 
Perron-Frobenius-Ruelle operator) . Continuing J5| , it is shown in [SJ that a gen- 
eral family of Hilbert-space constructions, which includes the context of wavelets, 
may be understood from the so-called Perron-Frobenius eigenspace of R. This 
eigenspace is studied further in the book JU] by Jorgensen and O. Bratteli; see 
also |55|. It is shown in |54| that a scaling identity (alias refinement equation) 

ip(x) = VN^2a k <p(Nx-k) (2.2) 

may be viewed generally and abstractly. Variations of (|2.2|) open up for a variety 
of new applications which we outline briefly below. 

3. Motivating examples, nonlinearity 

In this paper we outline how basic operator-theoretic notions from wavelet theory 
(such as the succesful multiresolution construction) may be adapted to certain 
state spaces in nonlinear dynamics. We will survey and extend our recent papers 
on multiresolution analysis of state spaces in symbolic shift dynamics X(A), and 
on a class of substitution systems X(r) which includes the Julia sets; see, e.g., our 
papers |421 1431 1441 143] . Our analysis of these systems X starts with consideration 
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of Hilbert spaces of functions on X. But already this point demands new tools. 
So the first step in our work amounts to building the right Hilbert spaces. That 
part relies on tools from both operator algebras, and from probability theory (e.g., 
path-space measures and martingales). 

First we must identify the appropriate measures on X, and secondly, we 
must build Hilbert spaces on associated extension systems Xoc , called generalized 
solenoids. 

The appropriate measures fi on X will be constructed using David Ruelle's 
thcrmodynamical formalism |71j : We will select our distinguished measures /i on 
X to have minimal free energy relative to a certain function W on X. The rela- 
tionship between the measures fi on X and a class of induced measures on the 
extensions X x is based on a random- walk model which we developed in |44 | , such 
that the relevant measures on X^ are constructed as path-space measures on paths 
which originate on X . The transition on paths is governed in turn by a given and 
prescribed function W on X. 

In the special case of traditional wavelets, W is the absolute square of a so- 
called low-pass wavelet filter. In these special cases, the wavelet-filter functions 
represent a certain system of functions on the circle T = R/Z which we outline 
below; see also our paper [45] . Even our analysis in Jjj] includes wavelet bases 
on affine Cantor sets in M. d . Specializing our Julia-set cases X = X(r) to the 
standard wavelet systems, we get X — T, and the familiar approach to wavelet 
systems becomes a special case. 

According to its original definition, a wavelet is a function ip G L 2 (R) that 
generates an orthonormal basis under the action of two unitary operators: the 
dilation and the translation. In particular, 



Uf(x) = _/(£), Tf(x) = f(x-1), feL 2 (R),xeR. (3.1) 



One of the effective ways of getting concrete wavelets is to first look for 
a multiresolution analysis (MR A), that is to first identify a suitable telescoping 
nest of subspaces (V ra ) ne z of L 2 (R) that have trivial intersection and dense union, 
V n = UV n -i, and Vo contains a scaling function (p such that the translates of ip 
form an orthonormal basis for Vo- 

The scaling function will necessarily satisfy an equation of the form 



{U j T k t/j I j,k e Z} 
must be an orthonormal basis for L 2 (R), where 




(3.2) 



called the scaling equation (see also 112.211 1. 

To construct wavelets, one has to choose a low-pass filter 




(3.3) 
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and from this obtain the scaling function tp as the fixed point of a certain cascade 
operator. The wavelet tp (mother function) is constructed from tp (a father function, 
or scaling function) with the aid of two closed subspaces Vo , and Wq . The scaling 
function tp is obtained as the solution to a scaling equation (|3.2fl . and its integral 
translates generates Vo (the initial resolution), and tp with its integral translates 
generating Wo := ViOVo, where Vi is the next finer resolution subspace containing 
the initial resolution space Vo, and obtained from a one-step refinement of Vo- 

Thus the main part of the construction involves a clever choice for the low- 
pass filter mo that gives rise to nice orthogonal wavelets. 



3.1. MRAs in geometry and operator theory 

Let X be a compact Hausdorff space, and let r : X — > X be a finite-to-one contin- 
uous endomorphism, mapping X onto itself. As an example, r = r(z) could be a 
rational mapping of the Riemann sphere, and X could be the corresponding Julia 
set; or X could be the state space of a subshift associated to a 0-1 matrix. 

Due to work of Brolin [35] and Ruelle 1% , it is known that for these examples, 
X carries a unique maximal entropy, or minimal free-energy measure fi also called 
strongly r- invariant; see Lemma f4 1 1 . 31 f |TT |) . and equations i|4.2U|l and I|4.27[l - il4.28|) 
below. For each point the measure \i distributes the "mass" equally on the 

finite number of solutions y to r(y) = x. 

We show that this structure in fact admits a rich family of wavelet bases. 
Now this will be on a Hilbert space which corresponds to L 2 (R) in the familiar 
case of multiresolution wavelets. These are the wavelets corresponding to scaling 
x to Nx, by a fixed integer N, N > 2. In that case, X = T, the circle in C, and 
r(z) — z N . So even in the "classical" case, there is a "unitary dilation" from L 2 (X) 
to L 2 (R) in which the Haar measure on T "turns into" the Lebesgue measure on 
R. 

Our work so far, on examples, shows that this viewpoint holds promise for 
understanding the harmonic analysis of Julia sets, and of other iteration systems. 
In these cases, the analogue of L 2 (R) involves quasi- invariant measures /ioo on a 
space Xcyo , built from A" in a way that follows closely the analogue of passing from 
T to R in wavelet theory. But ^rx, is now a path-space measure. The translations 
by the integers Z on L 2 (R) are multiplications in the Fourier dual. For the Julia 
examples, this corresponds to the coordinate multiplication, i.e., multiplication 
by z on L 2 (X,n), a normal operator. We get a corresponding covariant system 
on Xoo, where multiplication by f(z) is unitarily equivalent to multiplication by 
f(r(z)). But this is now on the Hilbert space L 2 (Aoo), and defined from the path- 
space measure /x^. Hence all the issues that are addressed since the mid-1980's 
for L 2 (R)-wavelets have analogues in the wider context, and they appear to reveal 
interesting spectral theory for Julia sets. 
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3.2. Spectrum and geometry: wavelets, tight frames, and Hilbert spaces on Julia 
sets 

Background. Problems in dynamics are in general irregular and chaotic, lacking 
the internal structure and homogeneity which is typically related to group actions. 
Lacking are the structures that lie at the root of harmonic analysis, i.e., the torus 
or R d . 

Chaotic attractors can be viewed simply as sets, and they are often lacking 
the structure of manifolds, or groups, or even homogeneous manifolds. Therefore, 
the "natural" generalization, and question to ask, then depends on the point of 
view, or on the application at hand. 

Hence, many applied problems are typically not confined to groups; examples 
are turbulence, iterative dynamical systems, and irregular sampling. But there are 
many more. Yet harmonic analysis and more traditional MRA theory begins with 
Mr, or with some such group context. 

The geometries arising in wider contexts of applied mathematics might be 
attractors in discrete dynamics, in substitution tiling problems, or in complex 
substitution schemes, of which the Julia sets are the simplest. These geometries 
are not tied to groups at all. And yet they are very algorithmic, and they invite 
spectral-theoretic computations. 

Julia sets are prototypical examples of chaotic attractors for iterative discrete 
dynamical systems; i.e., for iterated substitution of a fixed rational function r{z), 
for z in the Riemann sphere. So the Julia sets serve beautifully as a testing ground 
for discrete algorithms, and for analysis on attractors. In our papers 01], 02] we 
show that multiscale/MRA analysis adapts itself well to discrete iterative systems, 
such as Julia sets, and state space for subshifts in symbolic dynamics. But so far 
there is very little computational harmonic analysis outside the more traditional 
context of R d . 

Definition 3.2.1. Let S be the Riemann sphere, and let r: S — > S be a rational 
map of degree greater than one. Let r n be the n'th iteration of r, i.e., the n'th 
iterated substitution of r into itself. The Fatou set F(r) of r is the largest open 
set U in S such that the sequence r™, restricted to U, is a normal family (in the 
sense of Montel). The Julia set X(r) is the complement of F(r). 

Moreover, X(r) is known to behave (roughly) like an attractor for the discrete 
dynamical system r n . But it is more complicated: the Julia set X(r) is the locus 
of expanding and chaotic behavior; e.g., X(r) is equal to the closure of the set of 
repelling periodic points for the dynamical system r n in S. 

In addition, X = X(r) is the minimal closed set X, such that \X\ > 2, and 
X is invariant under the branches of r , the inverses of r. 

While our prior research on Julia sets is only in the initial stages, there is 
already some progress. We are especially pleased that Hilbert spaces of martingales 
have offered the essential framework for analysis on Julia sets. The reason is that 
martingale tools are more adapted to irregular structures than are classical Fourier 
methods. 
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To initiate a harmonic analysis outside the classical context of the tori and 
of R d , it is thus natural to begin with the Julia sets. Furthermore, the Julia sets 
are already of active current research interest in geometry and dynamical systems. 

But so far, there are only sporadic attempts at a harmonic analysis for Julia 
sets, let alone the more general geometries outside the classical context of groups. 
Our prior work should prepare us well for this new project. Now, for Julia-set 
geometries, we must begin with the Hilbert space. And even that demands new 
ideas and new tools. 

To accommodate wavelet solutions in the wider context, we built new Hilbert 
spaces with the use of tools from martingale theory. This is tailored to applications 
we have in mind in geometry and dynamical systems involving Julia sets, fractals, 
subshifts, or even more general discrete dynamical systems. 

The construction of solutions requires a combination of tools that are some- 
what non-traditional in the subject. 

A number of wavelet constructions in pure and applied mathematics have a 
common operator-theoretic underpinning. It may be illustrated with the following 
operator system: Tt some Hilbert space; U: H — > H a unitary operator; V: T — > 
U {TL) a unitary representation. Here T is the 1-torus, T = R/Z, and U {TL) is the 
group of unitary operators in TL. 

The operator system satisfies the identity 

Vizy 1 UV(z) = zU, zel (3.4) 

Definition 3.2.2. We say that U is homogeneous of degree one with respect to the 
scaling automorphisms defined by { V (z) | z G T } if l|3.4|l holds. In addition, U 
must have the spectral type defined by Haar measure. We say that U is homoge- 
neous of degree one if it satisfies <|3.4|) for some representation V(z) of T. 

In the case of the standard wavelets with scale N, TL — L 2 (M), where in 
this case the two operators U and T are the usual iV-adic scaling and integral 
translation operators; see 1)3. 1|) above giving U for the case N = 2. 

If V C TL is a resolution subspace, then V C C/ _1 V . Set Wo ■= f7 _1 Vo9Vo, 
Qa := the projection onto Wo, Qk — U~ QqU , k £ Z, and 

oo 

V(z)= zk Qk, zeT. (3.5) 

k— — oo 

Then it is easy to check that the pair (U, V (z)) satisfies the commutation identity 
HI . 

Let be a Parseval frame in Wo- Then 

{ C/ fe T„i/> J \ k,neZ} (3.6) 

is a Parseval frame for TL. (Recall, a Parseval frame is also called a normalized 
tight frame.) 

Turning the picture around, we may start with a system (U, V (z)) which 
satisfies (|3.4f) . and then reconstruct wavelet (or frame) bases in the form l|3.6|l . 
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To do this in the abstract, we must first understand the multiplicity function 
calculated for { T n | n € Z } when restricted to the two subspaces Vo and Wo- But 
such a multiplicity function may be be defined for any representation of an abelian 
C*-algebra acting on TL which commutes with some abstract scaling automorphism 
(V(z)). 

This technique can be used for the Hilbert spaces associated to Julia sets, to 
construct wavelet (and frame) bases in this context. 

Wavelet filters in nonlinear models. These systems are studied both in the the- 
ory of symbolic dynamics, and in C*-algebra theory, i.e., for the Cuntz-Krieger 
algebras. 

It is known |25l I71| that, for these examples, X carries a unique maximal 
entropy, or minimal free-energy measure fi. The most general case when (X, r) 
admits a strongly r-invariant measure is not well understood. 

The intuitive property of such measures fi is this: For each point x in X, /i 
distributes the "mass" equally on the finite number of solutions y to r(y) = x. Then 
\x is obtained by an iteration of this procedure, i.e., by considering successively the 
finite sets of solutions y to r n (y) = x, for all n — 1,2,...; taking an average, and 
then limit. While the procedure seems natural enough, the structure of the convex 
set K(X, r) of all the r-invariant measures is worthy of further study. For the case 
when X = T the circle, \i is unique, and it is the Haar measure on T. 

The invariant measures are of interest in several areas: in representation the- 
ory, in dynamics, in operator theory, and in C* -algebra theory. Recent work of 
Dutkay-Jorgensen |48l 1441 l4"2~] focuses on special convex subsets of K(X,r) and 
their extreme points. This work in turn is motivated by our desire to construct 
wavelet bases on Julia sets. 

Our work on K(X,r) extends what is know in traditional wavelet analysis. 
While linear wavelet analysis is based on Lebesgue measure on R d , and Haar 
measure on the d-torus T d , the nonlinear theory is different, and it begins with 
results from geometric measure theory: We rely on results from |25| . |5()| . 
and 61 . And our adaptation of spectral theory in turn is motivated by work by 
Baggett et al. in [3[7|. 

Related work reported in |filj further ties in with exciting advances on graph 
C*-algebras pioneered by Paul Muhly, Iain Raeburn and many more in the C*- 
algebra community; see |HI , [IDl , [131, EH], EU, EU, ESI EH], IHSI , E3 • 

We now prepare the tools in outline, function systems and self-similar oper- 
ators which will be used in Subsection 14. II in the construction of specific Hilbert 
bases in the context of Julia sets. 

Specifically, we outline a construction of a multiresolution/wavelet analysis 
for the Julia set X of a given rational function r(z) of one complex variable, i.e., 
X = Julia(r). Such an analysis could likely be accomplished with some kind of 
wavelet representation induced by a normal operator T which is unitarily equiv- 
alent to a function of itself. Specifically, T is unitarily equivalent to r(T), i.e., 
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UT = r(T)U, for some unitary operator U in a Hilbert space Tt. (3.7) 

Even the existence of these representations would seem to be new and significant 
in operator theory. 

There are consequences and applications of such a construction: First we 
will get the existence of a finite set of some special generating functions nii on 
the Julia set X that may be useful for other purposes, and which are rather 
difficult to construct by alternative means. The relations we require for a system 
too, . . . , tojv_i of functions on X(r) are as follows: 

Jf m *(y) m ](y) h (y) — &ijh{x) for a.e.x € X(r), (3.8) 

yex(r) 

r(y)=x 

where h is a Perron-Frobenius eigenfunction for a Ruelle operator R defined on 
L°°(X(r)). Secondly, it follows that for each Julia set X, there is an infinite- 
dimensional group which acts transitively on these systems of functions. The gen- 
erating functions on X we are constructing are analogous to the more familiar 
functions on the circle T which define systems of filters in wavelet analysis. In 
fact, the familiar construction of a wavelet basis in H — L 2 (R) is a special case of 
our analysis. 

In the standard wavelet case, the rational function r is just a monomial, i.e., 
r{z) — z N where N is the number of frequency bands, and X is the circle T in the 
complex plane. The simplest choice of the functions too, . . . , m^-i and h in this 
case is rrik(z) = z k , z E T, < k < N — 1, and h(z) = 1, z € T. This represents 
previous work by many researchers, and also joint work between Jorgensen and 
Bratteli [T5| . 

Many applied problems are typically not confined to groups; examples are 
turbulence, iterative dynamical systems, and irregular sampling. But there are 
many more. 

Even though attractors such as Julia sets are highly nonlinear, there are 
adaptations of the geometric tools from the analysis on R to the nonlinear setting. 

The adaptation of traditional wavelet tools to nonlinear settings begins with 
an essential step: the construction of an appropriate Hilbert space; hence the mar- 
tingales. 

Jorgensen, Bratteli, and Dutkay have developed a representation-theoretic 
duality for wavelet filters and the associated wavelets, or tight wavelet frames. 
It is based on representations of the Cuntz algebra On (where N is the scale 
number); see [T2|, @H, EH- ^ s a by-product of this approach they get infinite 
families of incquivalcnt representations of Oat which are of independent interest 
in operator-algebra theory. 

The mathematical community has long known of operator-theoretic power in 
Fourier-analytic issues, e.g., Rieffel's incompleteness theorem jHS] for Gabor sys- 
tems violating the Nyquist condition, cf. |34|. We now feel that we can address a 
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number ol significant problems combining operator-theoretic and harmonic analy- 
sis methods, as well as broadening and deepening the mathematical underpinnings 
of our subject. 

3.3. Multiresolution analysis (MRA) 

One reason various geometric notions of multiresolution in Hilbert space (starting 
with S. Mallat) have proved useful in computational mathematics is that these 
resolutions are modeled on the fundamental concept of dyadic representation of 
the real numbers; or more generally on some variant of the classical positional 
representation of real numbers. In a dyadic representation of real numbers, the 
shift corresponds to multiplication by 2 _1 . 

The analogue of that in Hilbert space is a unitary operator U which scales 
functions by 2 _1 . If the Hilbert space Ti is L 2 (M.), then a resolution is a closed 
subspace Vo in Ti which is invariant under U with the further property that the 
restriction of U to Vo is a shift operator. Positive and negative powers of U then 
scale Ti into a nested resolution system indexed by the integers; and this lets us 
solve algorithmic problems by simply following the standard rules from number 
theory. 

However, there is no reason at all that use of this philosophy should be 
restricted to function spaces on the line, or on Euclidean space. Other geometric 
structures from discrete dynamics admit resolutions as well. 

Most MRA, Fourier multiresolution analysis (FMRA) and generalized mul- 
tiresolution analysis (GMRA) wavelet constructions in pure and applied math- 
ematics have a common operator-theoretic underpinning. Consider the operator 
system (Ti, U, V (z)) from [|3.4JI above. 

Proposition 3.3.1. Once the system (Ti, U, V) is given as in (|3.4|) . a scale of closed 
subspaces V n (called resolution subspaces) may be derived from the spectral sub- 
spaces of the representation T 3 z i— > V (z), i.e., subspaces (V n ) ne z such that 
V n C Vn+lj 

= {0}, U^™ * s dense in Ti, 

and UV n C V n -i- 

Conversely, if these spaces (V n ) are given, then V(z) defined by equation 
(|3.5(l can be shown to satisfy (|3.4|) if 

Qo = the orthogonal projection onto C/ _1 Vo Vo, 

and Q n = U- n Q U n , neZ. 

As a result we note the following criterion. 

Proposition 3.3.2. A given unitary operator U in a Hilbert space Ti is part of some 
multiresolution system (V„) for Ti if and only if U is homogeneous of degree one 
with respect to some representation V(z) of the one-torus T. 
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In application the spaces V n represent a grading of the entire Hilbert space 
H, and we say that U scales between the different levels. In many cases, this 
operator-theoretic framework serves to represent structures which are similar up 
to scale, structures that arise for example in the study of fractional Brownian 
motion (FBM). See, e.g., which offers a Hilbert-space formulation of FBM 
based on a white-noise stochastic integration. 

Moreover, this can be done for the Hilbert space Ti — L 2 (^oo Moo) associated 
with a Julia set X. As a result, we get wavelet bases and Parseval wavelets in this 
context. 

This is relevant for the two mentioned classes of examples, Julia sets X(r), 
and state spaces X(A) of substitutions or of subshifts. Recall that the translations 
by the integers Z on L 2 (M) are multiplication in the Fourier dual. For the Julia 
sets, this corresponds to coordinate multiplication, i.e., multiplication by z on the 
Hilbert space L 2 (A, /i), consequently a normal operator. The construction will 
involve Markov processes and martingales. 

Thus, many of the concepts related to multiresolutions in L 2 (R) have ana- 
logues in a more general context; moreover, they seem to exhibit interesting con- 
nections to the geometry of Julia sets. 

Pyramid algorithms and geometry. Several of the co-authors' recent projects in- 
volve some variant or other of the Ruelle transfer operator R, also called the 
Ruelle-Perron-Frobenius operator. In each application, it arises in a wavelet-like 
setting. But the settings are different from the familiar L 2 (M)-wavelet setup: one 
is for affinc fractals, and the other for Julia sets generated from iteration dynamics 
of rational functions r(z). 

Thus, there are two general themes in this work. In rough outline, they are 
like this: (1) In the paper Dutkay and Jorgensen construct a new class of 
wavelets on extended fractals based on Cantor-like iterated function systems, e.g., 
the middle-third Cantor set. Our Hilbert space in each of our constructions is 
separable. For the fractals, it is built on Hausdorff measure of the same (fractal) 
dimension as the IFS-Cantor set in question. In the paper |45j. we further intro- 
duce an associated Ruelle operator R, and it plays a central role. For our wavelet 
construction, there is, as is to be expected, a natural (equilibrium) measure v 
which satisfies vR = i.e., a left Perron-Frobenius eigenvector. It corresponds 
to the Dirac point-measure on the frequency variable to = (i.e., low-pass) in 
the standard L 2 (R)-wavelet setting. It turns out that our measures v for the IFS- 
Cantor sets are given by infinite Riesz products, and they are typically singular 
and have support = the whole circle T. This part of our research is related to 
recent, completely independent, work by Benedetto et al. on Riesz products. 

(2) A second related research direction is a joint project in progress with 
Ola Bratteli, where we build wavelets on Hilbert spaces induced by Julia sets 
of rational functions r(z) of one complex variable. Let r be a rational function 
of degree at least 2, and let X(r) be the corresponding Julia set. Then there is 
a family of Ruelle operators indexed by certain potentials, and a corresponding 
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family of measures v on X(r) which satisfy vR = v, again with the measures 
v playing a role somewhat analogous to the Dirac point-mass on the frequency 
variable u = 0, for the familiar L 2 (R)-MRA wavelet construction. 

Independently of Dutkay-Jorgensen , John Benedetto and his co-authors 
Erica Bernstein and Ioannis Konstantinidis |12| have developed a new Fourier/in- 
finite-product approach to the very same singular measures that arise in the study 
in 132]. The motivations and applications are different. In 05] the issue is wavelets 
on fractals, and in |12| . it is time- frequency duality. 

In this second class of problems, the measures v are far less well understood; 
and yet, if known, they would yield valuable insight into the spectral theory of 
r-iteration systems and the corresponding Julia sets X(r). 

3.4. Julia sets from complex dynamics 

Wavelet-like constructions are just under the surface in the following brief sketch 
of problems. The projects all have to do with the kind of dynamical systems X 
already mentioned; and f2 is one of the standard probability spaces. The space X 
could be the circle, or a torus, or a solenoid, or a fractal, or it could be a Julia set 
defined from a rational map in one complex variable. More generally, consider maps 
r of X onto X which generate some kind of dynamics. The discussion is restricted 
to the case when for each x in X, the pre-image r~ 1 ({x}) is assumed finite. Using 
then a construction of Kolmogorov, it can be seen that these systems X admit 
useful measures which are obtained by a limit construction and by averaging over 
the finite sets r _1 ({a;}); see, e.g., Jorgensen's new monograph and Proposition 
I4.3.9l for more details. We list three related problems. 

(1) Operator algebras: A more systematic approach, along these lines, to 
crossed products with endomorphisms. For example, trying to capture the frame- 
work of endomorphism-crossed products introduced earlier by Bost and Connes 
|16| . (2) Dynamics: Generalization to the study of onto-maps, when the number 
and nature of branches of the inverse is more complicated; i.e., where it varies 
over X, and where there might be overlaps. (3) Geometry: Use of more detailed 
geometry and topology of Julia sets, in the study of the representations that come 
from multiplicity considerations. The Dutkay-Jorgensen paper 05] on fractals is 
actually a special case of Julia-set analysis. (One of the Julia sets X from the 
above is a Cantor set in the usual sense, e.g., the Julia set of r(z) = 2z — 1/z; or 
r{z) = z 2 - 3.) 

A main theme will be Hilbert spaces built on Julia sets X (r) built in turn 
on rational functions r (z) in one complex variable. In the case when r is not a 
monomial, say of degree N > 2, we construct N functions m,, i = 0, 1, . . . , iV — 1, 
on X (r) which satisfy a set of quadratic conditions 13.8fl analogous to the axioms 
that define wavelet filters. We have done this, so far, using representation theory, 
in special cases. But our results are not yet in a form where they can be used in 
computations. By using results on the Hausdorff dimension of Julia sets, it would 
also be interesting to identify inequivalent representations, both in the case when 
the Julia set X is fixed but the functions vary, and in the case when X varies. 
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Moreover, we get an infinite-dimensional "loop group" G — G(X, N) acting 
transitively on these sets of functions m,. The group G consists of measurable 
functions from X into the group Ujv of all unitary TV by N complex matrices. 
In particular, our method yields information about this group. Since X can be 
geometrically complicated, it is not at all clear how such matrix functions X — > 
Uat(C) might be constructed directly. 

The group G consists of measurable functions A: X(r) — > Ujv(C), and the 
action of A on an (TOi)-system is as follows: (mi) i— > (m\ ), where 

m\ A) {x) = J2A id ( r ( x )) m j( x )> x€X(r). (3.9) 

3 

Even when N = 2, the simplest non- monomial examples include the Julia sets 
of r(z) = z ! + c, and the two functions mo and mi on X(c) = Julia(z 2 + c) are not 
readily available by elementary means. The use of operator algebras, representation 
theory seems natural for this circle of problems. 

4. Main results 

4.1. Multiresolution/wavelet analysis on Julia sets 

We attempt to use equilibrium measures (Brolin [35] , Ruelle [71] , Mauldin-Urbanski 
61 ) from discrete dynamical systems to specify generalized low-pass conditions in 
the Julia-set theory. This would seem to be natural, but there are essential diffi- 
culties to overcome. While Ruelle's thermodynamical formalism was developed for 
transfer operators where the weight function W is strictly positive, the applica- 
tions we have in mind dictate careful attention to the case when W is not assumed 
strictly positive. The reason is that for our generalized multiresolution analysis on 
Julia sets, W will be the absolute square of the low-pass filter. 

We begin with one example and a lemma. The example realizes a certain 
class of state spaces from symbolic dynamics; and the lemma spells out a class 
of invariant measures which will be needed in our Hilbert-space constructions 
further into the section. The measures from part JuJ in Lemma 14.1.31 below are 
often called strongly invariant measures, or equilibrium measures from Ruelle's 
thermodynamical formalism, see |71| . 

Example 4.1.1. Let N <E Z+, N > 2 and let A = (ay)^- =1 be an TV by TV matrix 
with all a,ij e {0, 1}. Set 

X(A) := { (xA G . . . , N} | A( Xi , x i+1 ) = 1 } 

N 

and let r — be the restriction of the shift to X(A), i.e., 

rA(xi,x 2 ,---) = {x 2 ,x 3 , . . .), x = (xx,x 2 , ■ ■ ■) € X(A). 
Lemma 4.1.2. Let A be as above. Then 

#r^(x) = #{ye {!,..., N}\ A{y, x x ) = 1 }. 
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It follows that ta ■ X(A) — > X(A) is onto iff A is irreducible, i.e., iff for all 
j G Z^v, there exists an i £ Zjv such that A(i,j) = 1. 

Suppose in addition that A is aperiodic, i.e., there exists p G Z + such that 
A p > on Zjy x Zjv ■ We have the following lemma. 

Lemma 4.1.3. (D. Ruelle |72l IS]) Let A be irreducible and aperiodic and let (f> G 
C(X(A)) be given. Assume that (f> is a Lipschitz function. 

(i) Set 

(R<t,f)(x)= ]T e^f(y), forfeC(X(A)). 

r A (y)=x 

Then there exist Xq > 0. 

A = sup{ |A| | A G spec(i?0) }, 
h G C(X(A)) strictly positive and v a Borel measure on X(A) such that 

R<f>h = Xoh, 

and v{h) = 1. The data are unique. 

(ii) In particular, setting 

1 



(R f)(x) = : J2 

#r A (x) ^ 



r A (v)=x 

we may take Xq = 1, h = 1 and v —: /ia, where fiA is a probability measure 
on ^(^4) satisfying the strong invariance property 

fol^A^I - \, , E f(y) d »A(x), f£L°°(X(A). 

Our next main results, Theorems 14.1.41 and Theorem 14.2.11 are operator- 
theoretic, and they address the following question: Starting with a state space 
X for a dynamical system, what are the conditions needed for building a global 
Hilbert space H(X) and an associated multiscale structure? Our multiscales will 
be formulated in the context of Definitions 13 . 272\ for eg uivalent lv 13.3. lf> above; that 
is, we define our multiscale structure in terms of a natural Hilbert space H{X) 
and a certain unitary operator U on H{X) which implements the multiscales. 

After our presentation of the Hilbert-space context, we turn to spectral the- 
ory: In our next result, we revisit Baggett's dimension-consistency equation. The 
theorem of Baggett et al. [S] concerns a condition on a certain pair of multiplic- 
ity functions for the existence of a class of wavelet bases for L 2 (M. d ). In Theorem 
14.2.31 below, we formulate the corresponding multiplicity functions in the context 
of nonlinear dynamics, and we generalize the Baggett et al. theorem to the Hilbert 
spaces H(X) of our corresponding state spaces X. 

As preparation for Theorem 14.3.61 below, we introduce a family of _L 2 -mar- 
tingales, and we prove that the Hilbert spaces H(X) are L 2 -martingale Hilbert 
spaces. In this context, we prove the existence of an explicit wavelet transform 
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( Theorem 14 . 3 . 6|l for our nonlinear state spaces. But in the martingale context, we 
must work with dilated state spaces X^; and in Proposition ^. 3. 91 and Theorems 
15.8.11 15.8.21 15.8.91 and 15.8.121 we outline the role of our multiscale structure in 
connection with certain random- walk measures on X^. A second issue concerns 
extreme measures: in Theorem 15.8.81 we characterize the extreme points in the 
convex set of these random- walk measures. 

Theorem 4.1.4. Let A be a unital C* -algebra, a an endomorphism on A, fi a state 
on A and, mo £ A, such that 

H(m* a{f)m ) = fx{f), f e A. (4.1) 

Then there exists a Hilbert space H , a representation n of A on H , U a unitary 
on H , and a vector <p 6 A, with the following properties: 

Uir(f)U*=7r(a(f)), feA, (4.2) 



(vK./» = M/)> /^A ( 4 - 3 ) 

Uip = ir(a(l)m )(p (4.4) 

span{ U- n <ir(f)tp \ n>0JeA} = H. (4.5) 
Moreover, this is unique up to an intertwining isomorphism. 

We call (H, U, tt, tp) the covariant system associated to /i and mo- 

Corollary 4.1.5. Let X be a measure space, r: X — ► X a measurable, onto map 
and \i a probability measure on X such that 

I f ^ = I Wr^Ix) £ f(y)M*Y (4-6) 
Jx Jx#r {.x) r(yy=x 

Let h e L 1 (X), h>0 such that 

# r -i< x \ l m o(y)l 2 %) = H x )> xex. 

r{y)=x 

Then there exists {uniquely up to isomorphism) a Hilbert space H , a unitary U , a 
representation ir of L°°(X) and a vector tp € H such that 

UnifjU- 1 =7r(/or), / £ L°°(X), 

(tp\n(f)<p) = [ fhdfi, f&L°°(X), 
Jx 

Uip = ir(m Q )tp, 
span{ U- n n{f)cp \ n > 0, / € L°°{X) } = H. 
We call (H, U,tt, ip) the covariant system associated to nig and h. 
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4.2. The Baumslag Solitar group 

Our next theorem is motivated by the theory of representations of a certain group 
which is called the Baumslag-Solitar group, and which arises by the following 
simple semidirect product construction: the Baumslag-Solitar group BS(1, N) is 
the group with two generators u and t and one relation utu^ 1 — t N , where N is a 
fixed positive integer. Therefore the two operators U of ./V-adic dilation and T of 
integral translation operators on K give a representation of this group. But so do 
all familiar MRA-wavelet constructions. 

Representations of the Baumslag-Solitar group that admit wavelet bases can 
be constructed also on some other spaces, such as L 2 (K) © • • • © L 2 (R) (131 141| or 
some fractal spaces |4*5] . 

Hence wavelet theory fits naturally into this abstract setting of special rep- 
resentations of the group BS(1, N) realized on some Hilbert space. 

It is already known [3^1 that the representations that admit wavelets are 
faithful and weakly equivalent to the right regular representation of the group. 

Theorem 4.2.1. (i) Let H be a Hilbert space, S an isometry on H . Then there exist 
a Hilbert space H containing H and a unitary S on H such that 

S\ H = S, (4.7) 



|J S- n H = H. (4.8) 

n>0 

Moreover these are unique up to an intertwining isomorphism. 

(ii) If A is a C* -algebra, a is an endomorphism on A and ir is a representa- 
tion of A on H such that 

STr(g) = n(a(g))S, g e A, (4.9) 
then there exists a unique representation it on H such that 

H9)\ H = n(g), geA, (4.10) 



Sn(g) = n(a(g))S, g e A. (4.11) 

Corollary 4.2.2. Let X, r, and \x be as in Corollary \4-1.5\ Let I be a finite or 
countable set. Suppose H: X — > B(£ 2 (I)) has the property that H(x) > for 
almost every x € X, and H^ G L X (X) for alli,j S /. Let Mq: X — > B(£ 2 (I)) such 
that x <— > ||Mo(x)|| is essentially bounded. Assume in addition that 

77 , l u x E M*(y)H(y)M (y) = H(x), for a.e. x £ X. (4.12) 

r(y)=x 

Then there exists a Hilbert space K , a unitary operator U on K , a representation 
tt of L°°(X) on K , and a family of vectors (tp,) £ K , such that 

m(g)U- 1 =7r(<7° r ), g e L°°(X), 
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(iPi | *(f)<p s ) = f fHa dn, i,jel, fe L°°(X), 
Jx 

span{7r(/)^ | n > 0,/ G L°°(X), i £ 1} = K. 

These are unique up to an intertwining unitary isomorphism. (All functions are 
assumed weakly measurable in the sense that x i— > ( £ | F(x)i] ) is measurable for 
all^n£t 2 (I).) 

Suppose now that H is a Hilbert space with an isometry S on it and with a 
normal representation tt of L°°(X) on H that satisfies the covariance relation 

S7r(g) = 7r(gor)S, geL^(X). (4.13) 

Theorem 14.2.11 shows that there exists a Hilbert space H containing H, a 
unitary S on H and a representation tt of L°°(X) on H such that 

(V n := S~ n (H)) n form an increasing sequence of subspaces with dense union, 

SI h = S, 
tt\ h — TT, 
57r(g) = Tt(gor)S. 

Theorem 4.2.3. (i) V% = 5 l_1 (_ff) is invariant for the representation tt. The mul- 
tiplicity functions of the representation tt on V\, and on Vq — H , are related by 

m Vl (x)= TO ^(y), xeX. (4.14) 

r(y)=x 

(ii) // W := Vi 9 Vq = S~ l H G H, then 

m Vo (x) + m Wo (x) = ^2 m v„(y), xeX. (4.15) 

r(y)=x 

Proof. Note that S maps Vi to Vq, and the covariance relation implies that the 
representation tt on V\ is isomorphic to the representation tt t : g <— > n(go r) on Vq. 
Therefore we have to compute the multiplicity of the latter, which we denote by 

777.1/ . 

By the spectral theorem there exists a unitary isomorphism J: H(= Vq) — > 
L 2 (X, my , ji) , where, for a multiplicity function m: X — > {0, 1, . . . , oo}, we use 
the notation 

L 2 (X,m,ri :={f:X^ U xeX C m ^ \ f(x) G C m{x \ [ \\f(x)f d»(x) < oo}. 

Jx 

In addition J intertwines tt with the representation of L°°(X) by multiplication 
operators, i.e., 

(JvtG?) J-\i)){x) = g(x)f(x), g G L°°(X)J G L 2 (X, 777 y , /j), x G X 
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Remark 4.2.4. Here we are identifying H with L 2 (X, mv , m) y i a the spectral repre- 
sentation. We recall the details of this representation H 3 f i— > / S L 2 (X, my , A*)- 
Recall that any normal representation it £ Rep(L°°(X), H) is the orthogonal 

sum 

H=J2 [n(L°°(X))k], (4.16) 
kec 

where the set C of vectors k £ H is chosen such that 

11*11=1, 

• (k\ir(g)k) = / g(x)vk{x) 2 dfj,(x), for all fc £ C; 

J x 

• ( k' | vr(g)fc) =0, 5 £ L°°(X),k, k' £C,k^ orthogonality. 

The formula l|4.16|l is obtained by a use of Zorn's lemma. Here, v\ is the 
Radon-Nikodym derivative of {k | 7r(-)fc) with respect to fj,, and we use that n is 
assumed normal. 

For f £ H, set 

f=Y^<9k)k, g k £L™(X) 



kec 



and 



/= Y!" 9kVkeLl(X,£ 2 (C)). 



kec 



Then = / is the desired spectral transform, i.e., 

W is unitary, 
Ww(g) = M(g)W, 



and 



Indeed, we have 



l/»ll 2 = E MxHWf 



kec 



f \\f(x)\\ 2 d^(x) = f ^2\g k ( X )\ 2 v k (xfdn(x) = Y^ [ \gk\ 2 v 2 k dn 



kec 



H 



= E( fc i 7r (i^i 2 ) fc ) = En^)fcii 2 

kec kec kec 

It follows in particular that the multiplicity function m(x) = muix) is 

m(x) = #{k£C\v k {x)^0}. 

Setting 

Xi := { x £ X | m{x) > i }, i > 1, 

we see that 

~ E £ 2 (X,^)~L 2 (X,m, M ), 
and the isomorphism intertwines ir(g) with multiplication operators. 



\f\\l- 
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Returning to the proof of the theorem, we have to find the similar form for 
the representation tt t . Let 

m(x):= m v (y), x^X. (4.17) 

r(y)=x 

Define the following unitary isomorphism: 

L : L 2 (X,m Vo , g.) — > L 2 (X,m, fi), 

(L£)(x) = 1 = {€{y))r(y)=x- 

(Note that the dimensions of the vectors match because of (|4.17|) .) This operator 
L is unitary. For £ € L 2 (X, my 01 /i), we have 

ll^ll^(x,w ^) = / ll^)fW 

= / HC^II 2 ^). 

And L intertwines the representations. Indeed, for g £ L°°(X), 

L(por£)(x) = (g(r(y))£(y)) r ( y)=x = g(x)L(^)(x). 

Therefore, the multiplicity of the representation 7r r : g i— > 7r(g or) on Vo is m, and 
this proves (i). 

(ii) follows from (i). 

Conclusions. By definition, if fc € C, 

(fc|7r(ff)fc)= / g(x)v k (x) 2 dg,(x), and 

( I K r {g)k) = [ g(r(x))v k (x) 2 dg(x) = / g(x) — _\ V u fe (a;) 2 d/i(a;); 
Jx Jx w r \ x ) , s 



and so 

r(y)=x 

= £ #{fceC| Ufe (y) 2 >0} 

r(y)=x 
r(|/)=x 

Let C rn (x) := {k e C \ v k (x) ^0}. Then we showed that 

C m (x) = (J C m (y) 

2/£X,r(i/)=x 



Multiscale Theory and Wavelets in Nonlinear Dynamics 



21 



and that C m (y) n C m (y') = when y ^ y' and r{y) = r{y r ) = x. Setting H(x) = 
£ 2 (C m (x)), we have 

H{x) = f(C m (x)) = ]T e f(C m (y)) = ]T® H{y). □ 

r(y)=x r(y)=x 

4.3. Spectral decomposition of covariant representations: projective limits 

We give now a different representation of the construction of the covariant system 
associated to m and h. 

We work in either the category of measure spaces or topological spaces. 

Definition 4.3.1. Let r: X — ► X be onto, and assume that #r _1 (x) < oo for all 
x € X. We define the projective limit of the system: 

X^-X^X---< — X x (4.18) 

as 

Xoo := {x = (x ,xi, . . . ) | r(x n+1 ) = x n , for all n > 0}. 

Let 6 n : X^ — > X be the projection onto the n-th component: 

9 n (x ,xi, ...)= x n , {x ,xi, . . .) e -Xoo. 

Taking inverse images of sets in X through these projections, we obtain a sigma- 
algebra on X^ , or a topology on X^ . 

We have an induced mapping f : X x — > defined by 

r(x) = (r(a;o), ^o, £ij • • • ) ; and with inverse f _1 (i) = (x\, X2, . . . ). (4.19) 

so f is an automorphism, i.e., f of -1 = idx^ and r -1 o f = idx^ ■ 
Note that 

9 n o f = r o 6 n = n -i, 




For mo € define 

(flfl(s) = — L- ^ |m (2/)| 2 e(y), ^i'(I). (4.21) 
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We now resume the operator-theoretic theme of our paper, that of extending 
a system of operators in a given Hilbert space to a "dilated" system in an ambient 
or extended Hilbert space, the general idea being that the dilated system acquires 
a "better" spectral theory: for example, contractive operators dilate to unitary 
operators in the extended Hilbert space, and similarly, endomorphisms dilate to 
automorphisms. This, of course, is a general theme in operator theory; see, e.g., 
and |88| . But in our present setting, there is much more structure than the 
mere Hilbert-space geometry. We must adapt the underlying operator theory to 
the particular function spaces and measures at hand. The next theorem (Theorem 
I4.3.4f> is key to the understanding of the dilation step as it arises in our context 
of multiscale theory. It takes into consideration how we arrive at our function 
spaces by making dynamics out of the multiscale setting. Hence, the dilations we 
construct take place at three separate levels, as follows: 

• dynamical systems, 

(X,r,n) endomorphism — > (Xoo,f,(l), automorphism ; 

• Hilbert spaces, 

L 2 (X,hdfi) -> (R mo h = h) -> ^{X^, fi); 

• operators, 

S mo isometry — > U unitary (if mo is non-singular) , 

M(g) multiplication operator — > M^g). 

Definition 4.3.2. A function mo on a measure space is called singular if mo = on 
a set of positive measure. 

In general, the operators S mo on Hq = L {X,hdft), and U on L 2 (X oc ,/t), 
may be given only by abstract Hilbert-space axioms; but in our martingale repre- 
sentation, we get the following two concrete formulas: 

(Uf)(x) = m (x )f(r(x)), x € Xoo, f G L^X^ft). 

Definition 4.3.3. Let X, r, fi, and W be given as before. Recall r: X — » X is a 
finite-to-one map onto X, [i is a corresponding strongly invariant measure, and 
W : X — > [0, oo ) is a measurable weight function. We define the corresponding 
Ruelle operator R = Rw, and we let h be a chosen Perron-Frobenius-Ruelle 
eigenf unction. From these data, we define the following sequence of measures to n 
on X. (Our presentation of a wavelet transform in Theorem 14.3.61 will depend on 
this sequence.) For each n, the measure lu„ is defined by the following formula: 

= / R n {fh)dn, f S L°°(X). (4.22) 
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To build our wavelet transform, we first prove in Theorem 14.3.41 below that 
each of the measure families (u> n ) defines a unique VK-quasi-invariant measure fx on 
Xoo. In the theorem, we take W := |too| 2 , and we show that jx is quasi-invariant 
with respect to r with transformation Radon-Nikodym derivative W . 

Theorem 4.3.4. If h E L l (X), h > and Rh = h, then there exists a unique 
measure fx on Xoo such that 

fx o Q^ 1 = uj n , n > 0. 
We can identify functions on X with functions on X^ by 
f(xo,X!,...) = f{x ), /:X->C. 
Under this identification, 

■■ \m \ 2 . (4.23) 



dfx 

Theorem 4.3.5. Suppose mo is non-singular, i.e., it does not vanish on a set of 
positive measure. Define U on L 2 (X 00 ,/i) by 

Uf = m fof, f eL 2 (X 00 ,fx), 
n(9)f = 9f, geL 00 (X),feL 2 (X 00 ,fi), 
<p = l. 

Then (L 2 (Xoo, fx), U, it, ip) is the covariant system associated to too and h as in 
Corollary \4-1.5\ Moreover, if M g f = gf for g 6 L°°(X 00 , fx) and f £ L (Xoo, fx), 
then 

UMgU- 1 = M g0 f. 

The Hilbert space L 2 (X co , fx) admits a different representation as an ^-mar- 
tingale Hilbert space. Let 

H n :={/£LK« I / = £ ° Qn,£ € L 2 (X,uj n ) }. 

Then H n form an increasing sequence of closed subspaces which have dense union. 
We can identify the functions in H n with functions in L 2 (X, u> n ), by 

in(O=Z°0n, £,EL 2 {X,LU n ). 

The definition of ft makes i n an isomorphism between H n and L 2 (X, ui n ). 
Define 

H := | (Co, Si,-..) | £n € L 2 (X,u n ), R(£ n+ ih) = £ n h, sup^iT(|£ n | 2 /i) dii < ocj 
with the scalar product 

((£o,£i,---) I im,m, •••))= lim / R n (^ n Vnh)dLx. 

n ^°° Jx 



24 



Dutkay and Jorgensen 



Theorem 4.3.6. The map $: L 2 (Xoo,/i) — ► TL defined by 

Hf) = Vn/))n>0, 

where P n is the projection onto H n , is an isomorphism. The transform $ satisfies 
the following three conditions, and it is determined uniquely by them: 

$f7$ _:L (£„)n>o = (mo°r"^+i)«>o, 
$7r(ff)$ _1 (f n ) n > = (.gor n £„)„>o, 
^=(1,1,...). 

Theorem 4.3.7. There exists a unique isometry ^> : L^iX^, fi) — > H such that 

^ intertwines the two systems, i.e., 

vjr[/ = U^, ^Tr(g) = Tt(g)®, for g £ L°°(X, fi), ^ip = (p. 

Theorem 4.3.8. Let (mo,h) be a Perron-Ruelle-Frobenius pair with mo non-sin- 
gular. 

(i) For each operator A on L 2 (Xoo, /2) which commutes with U and ir, there exists 
a cocycle f, i.e., a bounded measurable function f: X^ — > C with f = / o f, 
fi-a.e., such that 

A = M f , (4.24) 
and, conversely each cocycle defines an operator in the commutant. 

(ii) For each measurable harmonic function ho: X — > C, i.e., R mo ho — ho, with 
\ho\ < eft 2 for some c > 0, there exists a unique cocycle f such that 

h = E (f)h, (4.25) 

where Eq denotes the conditional expectation from L°°(X !X> , p,) to L°°(X,pi), 
and conversely, for each cocycle the function ho defined by <|4.25[) is harmonic. 

(iii) The correspondence ho — > / in (JuJ is given by 

f = lim ^ o 6 n (4.26) 
where the limit is pointwise fi-a.e., and in L p (Xoo,/t) for all 1 < p < oo. 

The next result in this section concerns certain path-space measures, indexed 
by a base space X. Such a family is also called a process: it is a family of positive 
Radon measures P x , indexed by x in the base space X . Each P x is a measure on the 
probability space Q which is constructed from X by the usual Cartesian product, 
i.e., countably infinite Cartesian product of X with itself. The Borel structure on 
ft is generated by the usual cylinder subsets of ft. Given a weight function W 
on X, the existence of the measures P x comes from an application of a general 
principle of Kolmogorov. 

Let X be a metric space and r : X — > X an N to 1 map. Denote by : X — » 
X, k £ {1, . . . , N}, the branches of r, i.e., r(rfe(x)) = x for x £ X, the sets Tk(X) 
are disjoint and they cover X . 
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Let n be a measure on X with the property 

1 - 

k=l 

This can be rewritten as 

N 

I f(x)d,x(x) = ±-J2 I f(Tk(x))d»(x), (4.28) 
Jx fc=lV X 

which is equivalent also to the strong invariance property. 
Let W, h > be two functions on X such that 

N 

£ W(r k (x))h{r k (xj) = h(x), xeX. (4.29) 
fe=l 

Denote by f2 the multi-index set 

n ~Q N := JJ{1,...,JV}. 

N 

Also we define 

W {n) (x) := W(x)W(r(x)) ■ ■ ■ W^r™" 1 ^)), x e X. 

We will return to a more general variant of these path-space measures in 
Theorem l5.8.12l in the next section. This will be in the context of a system (X, r, V) 
where r: X — > X is a give finite to one endomorphism, and W: X — * [0, oo) is a 
measurable weight function. In this context the measures (P x ) are indexed by a 
base space X, and they are measures on the projective space X^ built on [X,r), 
see Definition ^. 3. II The family (P x ) is also called a process: each P x is a positive 
Radon measure on Xoo, and X^ may be thought of as a discrete path space. 

The Borel structure on O, and on X^, is generated by the usual cylinder 
subsets. Given a weight function W on X, the existence of the measures P x := 
comes from an application of a general principle of Kolmogorov. 

Proposition 4.3.9. For every x G X there exists a ■positive Radon measure P x on 
n such that, if f is a bounded measurable function on X which depends only on 
the first n coordinates ui,..., Lo n , then 

f f{uj)dP x (uj) 
J a 

= W {n \T^ n _ 1 ---T^{x))h{T Un T Un _ 1 ---T^ 1 {x))f{u0 1 ,...,UJ n ). (4.30) 

UJl ,. . . ,UJ n 
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5. Remarks on other applications 

5.1. Wavelet sets 

Investigate the existence and construction of wavelet sets and elementary wavelets 
and frames in the Julia-set theory, and the corresponding interpolation theory, and 
their relationships to generalized multiresolution analysis (GMRA). The unitary 
system approach to wavelets by Dai and Larson [321 dealt with systems that can 
be very irregular. And recent work by Olafsson et al. shows that wavelet sets can 
be essential to a basic dilation-translation wavelet theory even for a system where 
the set of dilation unitaries is not a group. 

5.2. The renormalization question 

When are there renormalizable iterates of r in arithmetic progression, i.e., when are 
there iteration periods n such that the system {r kn , k € N} may be renormalized, 
or rescaled to yield a new dynamical system of the same general shape as that of 
the original map r(z)7 

Since the scaling equation (|2.2I) from wavelet theory is a renormalization, 
our application of multiresolutions to the Julia sets X(r) of complex dynamics 
suggests a useful approach to renormalization in this context. The drawback of 
this approach is that it relies on a rather unwieldy Hilbert space built on X(r), 
or rather on the projective system X(r)oa built in turn over X(r) . So we are left 
with translating our Hilbert-space-theoretic normalization back to the direct and 
geometric algorithms on the complex plane. 

General and rigorous results from complex dynamics state that under addi- 
tional geometric hypotheses, renormalized dynamical systems range in a compact 
family. 

The use of geometric tools from Hilbert space seems promising for renormal- 
ization questions (see, e.g., [S3], JS], [23], and [73]) since notions of repetition up 
to similarity of form at infinitely many scales are common in the Hilbert-space 
approach to multiresolutions; see, e.g., [35]. And this self-similarity up to scale 
parallels a basic feature of renormalization questions in physics: for a variety of in- 
stances of dynamics in physics, and in telecommunication |64lll5ll2TT] . we encounter 
scaling laws of self-similarity; i.e., we observe that a phenomenon reproduces itself 
on different time and/or space scales. 

Self-similar processes are stochastic processes that are invariant in distribu- 
tion under suitable scaling of time and/or space (details below!) Fractional Brown- 
ian motion is perhaps the best known of these, and it is used in telecommunication 
and in stochastic integration. While the underlying idea behind this can be traced 
back to Kolmogorov, it is only recently, with the advent of wavelet methods, that 
its computational power has come more into focus, see e.g., [S3|. But at the same 
time, this connection to wavelet analysis is now bringing the operator-theoretic 
features in the subject to the fore. 

In statistics, we observe that fractional Brownian motion (a Gaussian pro- 
cess B(t) with E{B(t)} = 0, and covariance E{B(t)B(s)} given by a certain 
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h- fractional law) has the property that there is a number h such that, for all a, 
the two processes B(at) and a h B{t) have the same finite-dimensional distribu- 
tions. The scaling feature of fractional Brownian motion, and of its corresponding 
white-noise process, is used in telecommunication, and in stochastic integration; 
and this connection has been one of the new and more exciting domains of appli- 
cations of wavelet tools, both pure and applied (Donoho, Daubechies, Meyer, etc. 

1331321 EES). 

This is a new direction of pure mathematics which makes essential contact 
with problems that are not normally investigated in the context of harmonic anal- 
ysis. 

In our proofs, we take advantage of our Hilbert-space formulation of the 
notion of similarity up to scale, and the use of scales of closed subspaces V„ in 
a Hilbert space Ti. The unitary operator U which scales between the spaces may 
arise from a substitution in nonlinear dynamics, such as in the context of Julia 
sets [TT]; it may "scale" between the sigma- algebras in a martingale |3o1 1571 175] : 
it may be a scaling of large volumes of data (see, e.g., U); it may be the scaling 
operation in a fractional Brownian motion 64 59 ; it may be cell averages in finite 
elements [201; or ^ may be the traditional dyadic scaling operator in Tt = L 2 (R) 
in wavelet frame analysis |28l l2"§] . 

Once the system (7i, (V n ), U) is given, then it follows that there is a unitary 
representation V of T in Ti which defines a grading of operators on TL. Moreover 
U then defines a similarity up to scale precisely when U has grade one, see (|3.4|) 
for definition. We also refer to the papers 0J E3 for details and applications of 
this notion in operator theory. 

The general idea of assigning degrees to operators in Hilbert space has served 
as a powerful tool in other areas of mathematical physics (see, e.g., |3J, |22]; an( i 
[TB| on phase transition problems); operator theory [5J 0J |3H1 : an( A operator 
algebras; see, e.g., |5T1 1751 12"4"] . 

Our research, e.g., |43j . in fact already indicates that the nonlinear problems 
sketched above can be attacked with the use of our operator-theoretic framework. 

5.3. Wavelet representations of operator-algebraic crossed products 

The construction of wavelets requires a choice of a low-pass filter mo <G L°° (T) and 
of N— 1 high-pass filters mi G L°°(T) that satisfy certain orthogonality conditions. 
As shown in J7J, these choices can be put in one-to-one correspondence with a 
class of covariant representations of the Cuntz algebra On- In the case of an 
arbitrary iV-to-one dynamical system r: X — > X, the wavelet construction will 
again involve a careful choice of the filters, and a new class of representations is 
obtain. The question here is to classify these representations and see how they 
depend on the dynamical system r. A canonical choice of the filters could provide 
an invariant for the dynamical system. 

In a recent paper 07], the crossed-product of a C*-algebra by an endomor- 
phism is constructed using the transfer operator. The required covariant relations 
are satisfied by the operators that we introduced in 03] , |44|. 
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There are two spaces for these representations: one is the core space of the 
multiresolution Vq — L 2 (X, fi) . Here we have the abelian algebra of multiplication 
operators and the isometry S: / i— » mo for. Together they provide representations 
of the crossed product by an endomorphism. The second space is associated to a 
dilated measure fx on the solenoid of r, L 2 (X 00 ,/x) (see 021; HH)- The isometry 
S is dilated to a unitary S which preserves the covariance. In this case, we are 
dealing with representations of crossed products by automorphisms. 

Thus we have four objects interacting with one another: crossed products by 
endomorphisms and their representations, and crossed products by automorphisms 
and its representations on L 2 (X,ft). 

The representations come with a rich structure given by the multiresolution, 
scaling function and wavelets. Therefore their analysis and classification seems 
promising. 

5.4. KMS-states and equilibrium measures 

In [38], the KMS states on the Cuntz-Pimsner algebra for a certain one-parameter 
group of automorphisms are obtained by composing the unique equilibrium mea- 
sure on an abelian subalgebra (i.e., the measure which is invariant to the transfer 
operator) with the conditional expectation. These results were further generalized 
in [HE] for expansive maps and the associated groupoid algebras. 

As we have seen in |45] for the case when r{z) = z N , the covariant representa- 
tions associated to some low-pass filter mo are highly dependent on the equilibrium 
measure v. As outlined above, we consider here an expansive dynamical system 
built on a given finite-to-one mapping r: X — > X, and generalizing the familiar 
case of the winding mapping z — ► z N on T. In our r: X — > X context, we then 
study weight functions V : X — > [0, oo). If the zero-set of V is now assumed finite, 
we have so far generalized what is known in the z — > z N case. Our results take the 
form of a certain dichotomy for the admissible Perron-Frobenius-Ruelle measures, 
and they fit within the framework of |58| . Hence our construction also yields covari- 
ant representations and symmetry in the context of |58| . Our equilibrium measures 
in turn induce the kind of KMS-states studied in [55] . 

The results of [SHj and 0H| restrict the weight function (in our case repre- 
sented by the absolute square of the low-pass filter, |mo| 2 ) to being strictly pos- 
itive and Holder continuous. These restrictions are required because of the form 
of the Ruelle-Perron-Frobenius theorem used, which guarantees the existence and 
uniqueness of the equilibrium measure. 

However, the classical wavelet theory shows that many interesting examples 
require a "low-pass condition", m (l) = VN, which brings zeroes for |m | 2 . 

Our martingale approach is much less restrictive than the conditions of |58| : 
it allows zeroes and discontinuities. Coupled with a recent, more general form of 
Ruelle's theorem from 0§| , we hope to extend the results of |2H] and be able to find 
the KMS states in a more general case. Since the existence of zeroes can imply 
a multiplicity of equilibrium measures (see [41| . \W\ ) a new phenomenon might 
occur such as spontaneously breaking symmetry. 
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For each KMS state, one can construct the GNS representation. Interest- 
ing results are known j^S] for the subshifts of finite type and the Cuntz-Krieger 
algebra, when the construction yields type IIIa AFD factors. 

The natural question is what type of factors appear in the general case of 
r: X — > X. Again, a canonical choice for the low-pass filter can provide invariants 
for the dynamical system. 

Historically, von-Neumann-algebra techniques have served as powerful tools 
for analyzing discrete structures. Our work so far suggests that the iteration 
systems {X, r) will indeed induce interesting von-Neumann-algebra isomorphism 
classes. 

5.5. Dimension groups 

We propose to use our geometric algorithms 20, 21 for deciding order isomorphism 
of dimension groups for the purpose of deciding isomorphism questions which arise 
in our multiresolution analysis built on Julia sets. The project is to use the results 
in (20) to classify substitution tilings in the sense of [HE]; |67| . 

5.6. Equilibrium measures, harmonic functions for the transfer operator, and in- 
finite Riesz products 

In JHI, |41j . for the dynamical system z i-^ z N , we were able to compute the 
equilibrium measures, the fixed points of the transfer operator and its ergodic 
decomposition by analyzing the commutant of the covariant representations asso- 
ciated to the filter mo. In 03] j E2] we have shown that also in the general case 
there is a one-to-one correspondence between harmonic functions for the transfer 
operator and the commutant. Thus an explicit form of the covariant representa- 
tion can give explicit solutions for the eigenvalue problem R mo h — h, which are of 
interest in ergodic theory (see [§]). 

In some cases, the equilibrium measures are Riesz products (see these 
are examples of exotic measures which arise in harmonic analysis [121 1461 1561 |SJ 
1691 1141 1261 1271 1571 1621 ITU] . The wavelet-operator-theoretic approach may provide 
new insight into their properties. 

5.7. Non-uniqueness of the Ruelle Perron Frobenius data 

A substantial part of the current literature in the Ruelle-Perron-Frobenius oper- 
ator in its many guises is primarily about conditions for uniqueness of KMS; so 
that means no phase transition. This is ironic since Ruelle's pioneering work was 
motivated by phase-transition questions, i.e., non-uniqueness. 

In fact non-uniqueness is much more interesting in a variety of applications: 
That is when we must study weight functions W in the Ruelle operator R = Rw 
when the standard rather restrictive conditions on W are in force. The much cele- 
brated Ruelle theorem for the transition operator R gives existence and uniqueness 
of Perron-Frobenius data under these restrictive conditions. But both demands 
from physics, and from wavelets, suggest strong and promising research interest 
in relaxing the stringent restrictions that are typical in the literature, referring to 
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assumptions on W . There is a variety of exciting possibilities. They would give us 
existence, but not necessarily uniqueness in the conclusion of a new Ruelle-type 
theorem. 

5.8. Induced measures 

A basic tool in stochastic processes (from probability theory) involves a construc- 
tion on a "large" projective space X^, based on some marginal measure on some 
coordinate space X. Our projective limit space Xoo will be constructed from a 
finite branching process. 

Let A be a k x k matrix with entries in {0, 1}. Suppose every column in A 
contains an entry 1. 

Set X(A) := { (&)ieN € 1\ N {1, ...,k}\ A(&,&+i) = 1 } and 

rA(a,6,---) = (6,6,...) for£eX(A). 

Then ta is a subshift, and the pair {X{A),r a) satisfies our conditions. 

It is known that, for each A, as described, the corresponding system 
ta '■ X(A) — > X(A) has a unique strongly r^-invariant probability measure, pa, 
i.e., a probability measure on X(A) such that pa Ri = Pa, where R\ is defined 
as in H5.1fl below. 

We now turn to the connection between measures on X and an associated 
family of induced measures on Xoo , and we characterize those measures Xoc which 
are quasi-invariant with respect to the invertible mapping f, where X^ := { x = 
(xq,Xi, . ..) S IIn,, x I r ( x n+i) = x n }, f(x) = (r(x ),x ,x 1 , ...). 

In our extension of measures from X to X^, we must keep track of the 
transfer from one step to the next, and there is an operator which accomplishes 
this, Ruelle's transfer operator 

R w f(x)= 1 w (y)f(y)> feL\x,p). (5.i) 

In its original form it was introduced in |72j . but since that, it has found a variety 
of applications, see e.g., [H]. For use of the Ruelle operator in wavelet theory, we 
refer to [3] and 

The Hilbert spaces of functions on I m will be realized as a Hilbert spaces 
of martingales. This is consistent with our treatment of wavelet resolutions as 
martingales. This was first suggested in |51| in connection with wavelet analysis. 

In 02| ) we studied the following restrictive setup: we assumed that X carries 
a probability measure p. which is strongly r-invariant. By this we mean that 

/ fdp= I _\ f(y)Mx), feL°°(X). (5.2) 

Jx Jx#r (x) yeXr[y)=x 

If, for example X = R/Z, and r(x) = 2x mod 1, then the Haar measure on 
R/Z = Lebesgue measure on [0,1) is the unique strongly r-invariant measure on 
X. 
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Suppose the weight function V is bounded and measurable. Then define R = 
Rv, the Ruelle operator, by formula (|5.1|) . with V = W. 

Theorem 5.8.1. f|43|) Letr: X — > X and Xoo(r) be as described above, and suppose 
that X has a strongly r-invariant measure fj,. Let V be a non-negative, measurable 
function on X , and let Ry be the corresponding Ruelle operator. 

(i) There is a unique measure p on X^r) such that 

(a) Ao^o 1 «M(^/ l = ^p), 

(b) J x f dp o 9- 1 = J x R v (fh) dfi, n e N . 

(ii) The measure p on I m (r) satisfies 

d(p o f) 



dp 



= Vo6 , (5.3) 



Ryh = h. (5.4) 

If the function V : X — > [0, oo) is given, we define 

V {n \x) := V(z)V(r(a:)) • • • V(r n_1 (x)), 

and set dfj, n := V' n 'djUo- Our result states that the corresponding measure p on 
Xoo (r) is V-quasi- invariant if and only if 

d//o = (V d/io) o r . (5-5) 

Theorem 5.8.2. Let V: X — > [0, oo) be 2$ -measurable, and let fio be a measure on 
X satisfying the following fixed-point property 

d/io = (V duo) o r . (5-6) 

Then there exists a unique measure p on X^ (r) such that 

d(p o f ) 



dp 



Vo9 (5.7) 



p o 9 1 = ^o- 

Definition 5.8.3. Let V: X — > [0, oo) be bounded and 2J-measurable. We use the 
notation 

M V (X) :={/ie M{X) | d^ = {Vd^or- 1 }. 
For measures A on (Xoo(r), *Bco) we introduce 

M£tXoo(r)) := | A e M (Jfoo(r)) | A o f « A and = V o O J ■ (5.8) 

As in Definition 15.8.31 let X, r, and V be as before, i.e., r is a finite-to-one 
endomorphism of X, and V: X — > [0, oo) is a given weight function. In the next 
theorem, we establish a canonical bijection between the convex set M v (X) of 
measures on X with the set of y-quasi- invariant measures on X^ . which we call 
A#(*«,(r)), see £3. 

The results so far in this section may be summarized as follows: 



32 



Dutkay and Jorgensen 



Theorem 5.8.4. Let V be as in Definition \5.8.Sl For measures fj, on Xoo(r) and 
n € No, define 

c n {fi) -doe- 1 . 

Then Co is a bijective affine isomorphism of M^ i (X 00 (r)) onto M v (X) that pre- 
serves the total measure, i.e., Cq((1)(X) = /((^(r)) for all fi £ M^(X 00 (r)). 

Theorem 5.8.5. Let V: X — ► [0, oo) be continuous. Assume that there exist some 
measure v on (X, 03) and two numbers < a < b such that 

a < v(X) <b, anda< [ V (n) dv <b for all neN. (5.9) 
Jx 

Then there exists a measure p,Q on (X, 03) that satisfies 

dp Q = (VdfM)) o r _1 , 

and there exists a V -quasi-invariant measure fj on (X^Jr), 03oo). 

Theorem 5.8.6. Let (X, 03), and r: X — > X , be as described above. Suppose V : X — > 
[0, oo) is measurable, 

r(y)=x 

and that some probability measure vy on X satisfies 

vy ° Rv = vy. (5.10) 

Assume also that (X, 03) carries a strongly r-invariant probability measure p, such 
that 

p({x e X | V(x) > 0}) > 0. (5.11) 

Then 

(i) Ty(dp) = Ry(l)dp, for n £ N, where 1 denotes the constant function one. 

(ii) [Monotonicji!/] • ■ ■ < i?y +1 (l) < i?y(l) < ■ ■ ■ < 1, pointwise on X. 

(iii) The limit lim n _> R v (l) — hv exists, Ryhy — hy , and 

p{{xeX\ h v {x) > }) > 0. (5.12) 

(iv) The measure dfi V ^ = hy dp is a solution to the fixed-point problem 

1 v(p )-Mo • 

(v) The sequence dp^^ — hy dp defines a unique fS v ^ as in Theorem, 1 5. H.l\ 
and 

(vi) fJLn^{f) = J x Ry(fhy) dp for all bounded measurable functions f on X , and 
all neN. 

Finally, 

(vii) The measure fS v ^ on X^r) satisfying fS v " > o B^ 1 = //„ has total mass 

A°°(*oo(r))=p(M= / h v (x)dp(x). 

Jx 
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Definition 5.8.7. A measure (j,q € M 1 / (X) is called relatively ergodic with respect to 
(r, V) if the only non-negative, bounded 03-measurable functions / on X satisfying 

E no ( V f) = E no 00/ ° r i pointwise ° r -1 -a.e., 
are the functions which are constant /zo-a.e. 

Since we have a canonical bijection between the two compact convex sets of 
measures (X) and M^l 1 (X oa (r)), the natural question arises as to the extreme 
points. This is answered in our next theorem. We show that there are notions ergod- 
icity for each of the two spaces X and Xoc (r) which are equivalent to extremality 
in the corresponding compact convex sets of measures. 

Theorem 5.8.8. Let V: X — > [0, oo) be bounded and measurable. Let 
(i G Ml^X^r)), and /i :=/io O 1 e M?(X). 

The following assertions are eguivalent: 

(i) fj, is an extreme point of M^ i l (X ao (r)) ; 

(ii) V o 9 dfj, is ergodic with respect to r; 

(iii) no is an extreme point of MY(X); 

(iv) [1q is relatively ergodic with respect to (r, V) . 

We now turn to the next two theorems. These are counterparts of our di- 
mension-counting functions which we outlined above in connection with Theorem 
14.2.31 see especially Remark 14.2.41 They concern the correspondence between the 
two classes of measures, certain on X (see Theorem I5.8.8JI . and the associated 
induced measures on the projective space X^. Our proofs are technical and will 
not be included. (The reader is referred to 02 f° r further details.) Rather we 
only give a few suggestive hints: Below we outline certain combinatorial concepts 
and functions which are central for the arguments. Since they involve a counting 
principle, they have an intuitive flavor which we hope will offer some insight into 
our theorems. 

Let X be a non-empty set, 58 a sigma-algebra of subsets of X, and r : X — > X 
an onto, finite-to-one, and 93-measurable map. 

We will assume in addition that we can label measurably the branches of the 
inverse of r. By this, we mean that the following conditions are satisfied: 

The map c: X 3 x i— > #r _1 (a;) < oo is measurable . (5.13) 

We denote by A{ the set 

A, := {x e X | c(x) = #r" 1 (a;) > i}, i € N. 

Equation (|5.13() implies that the sets Ai are measurable. Also they form a decreas- 
ing sequence and, since the map is finite-to-one, 

oo 

X = \J(A l+1 \A l ). 

i=l 
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Then, we assume that there exist measurable maps r^: A4 — > X, i 6 {1,2,...} 
such that 

r~ 1 (x) = {7i(a!),... ) T ( ( x )(s)}, xeX, (5.14) 



r i (A i )e5B for alHe {1,2,...}. (5.15) 

Thus Ti(x), . . . , T c ( x )(a;) is a list without repetitions of the "roots" of x, r~ 1 (x). 
From l|5.14|l we obtain also that 

T i (A i )nr j (A j )=®, ifi^i, (5.16) 

and 

00 

'r i (A i )=X (5.17) 



u 

i=l 

In the sequel, we will use the following notation: for a function / : X 
we denote by / o n the function 



foTi(x) := 



f(n(x)) ifieA,, 
if x G X \ Ai 



Our Theorem 14 . 3 . 41 depends on the existence of some strongly invariant mea- 
sure /1 on X, when the system {X, r) is given. However, in the general measurable 
category, such a strongly invariant measure fi on X may in fact not exist; or if it 
does, it may not be available by computation. In fact, recent wavelet results (for 
frequency localized wavelets, see, e.g., [7] and [5]) suggest the need for theorems 
in the more general class of dynamical systems (X, r) . 

In the next theorem (Theorcm l5.8.9|) . we provide for each system, X, r, and 
V a substitute for the existence of strongly invariant measures. We show that there 
is a certain fixed-point property for measures on X which depends on V but not 
on the a priori existence of strongly r-invariant measures, and which instead is 
determined by a certain modular function A on X . This modular function in turn 
allows us to prove a dis-integration theorem for our ^-quasi-invariant measures 
// on Xoc . In Theorem 15.8.121 we give a formula for this dis-integration of a V- 
quasi-invariant measure fi on Xoo in the presence of a modular function A. Our 
dis-integration of /t is over a Markov process P x , for x in X, which depends on A, 
but otherwise is analogous to the process P x we used in Proposition ^. 3. 9l 

Theorem 5.8.9. Let (X, 23) and r: X — * X be as above. Let V: X — > [0, 00) be a 

bounded ^-measurable map. For a measure /io on {X, 58), the following assertions 
are equivalent. 

(i) The measure /iq has the fixed-point property 

Vfordfx = [ fd» , for all f G L°°(X, Mo ). (5.18) 
x J x 
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(ii) There exists a non-negative, 03 -measurable map A (depending on V and /io) 
such that 

A(j/) = 1, for fio-a.e. x e X, (5.19) 

r(y)=x 

and 

[ Vfdno= [ V A(y)/(y)dtt)(s), for all f e L°°(X,no). (5.20) 

« Jx , ■, 

r{y)=x 

Moreover, when the assertions are true, A is unique up to V dfiQ-measure zero. 

We recall the definitions: if 03 is a sigma-algebra on a set X and r : X — > X 
is a finite-to-one, onto and measurable map, we denote by Xoo the set 



^oc(r) := | {x ,x u ...) G [] X 



r(x n +i) = x ra , for all n G No 



We denote the projections by 8 n : X^r) 9 (xq,x\, . . .) i— > G X. The union 
of the pull-backs 0~ 1 (Q3) generates a sigma-algebra 03 oo. The map r extends to a 
measurable bijection f: Xoo(r) — ► Xoo(r) defined by 

f(x ,xi,...) = (r(x Q ),x ,x 1 ,...). 

Let V: X — * [0, oo) be a measurable, bounded function. We say that a mea- 
sure no on (X, 03) has the fixed-point property if 



Vford t i = fdfio, feL°°(X,iMi). 
x Jx 

We say that a measure /i on (X^r), 03^) is V -quasi-invariant if 

<i(/t o r) = V o O <^A- 
We recall the following result from |44j . 

Theorem 5.8.10. There exists a one-to-one correspondence between measures hq 
on X with the fixed-point property and V -quasi-invariant measures fi on X^r), 
given by 

fi Q = [i o Oq 1 , 

Proposition 5.8.11. Let (X, 03), r : X — > X and be as above, and let D: X — > [0, oo) 

be a measurable function with the property that 

E D (y) = l - ( 5 - 21 ) 

r(y)=x 

Denote by the product of compositions 

D [n) :=D-Dor Kor"" 1 , n G N, D (0) := 1. (5.22) 
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Then for each xq € X , there exists a Radon probability measure P Xa on Xo such 
that, if f is a bounded measurable function on Q xo which depends only on the first 
n coordinates x\, . . . , x n , then 

f(u)dP Xo {Lu)= D^\x n )f{ Xl ,...,x n ). (5.23) 

r n (x n )=x 

Theorem 5.8.12. Let (X, 03), r: X — > X and V: X — > [0, oo) be as above. Let /xo 
be a measure on (X, 03) with the fixed-point property (|5.18|l . Let A: X — > [0, 1] &e 
£/ie function associated to V and /iq as in Theorem \5.8.!A and let fi be the unique 
V -quasi-invariant measure on Xoo(r) with 

jlo9 1 = no, 

as in Theorem \5.8.1(A For A define the measures P Xo as in Proposition \5.8.11\ 
Then, for all bounded measurable functions f on X^r), 

/ fdp,= / / f(xo,u)dP xo (u)d(io(xo). (5.24) 
5.9. Hausdorff measure and wavelet bases 

In this problem we propose wavelet bases in the context of Hausdorff measure of 
fractional dimension between and 1. While our fractal wavelet theory has points 
of similarity that it shares with the standard case of Lebesgue measure on the line, 
there are also sharp contrasts. 

It is well known that the Hilbert spaces L 2 (R) has a rich family of orthonormal 
bases of the following form: ipj^(x) = 2^i 2 ip{2Px — k), j,k € Z, where -0 is a single 
function € L 2 (R). Take for example Haar's function ip(x) = x/(2x) — X/(2x — 1) 
where / = [0, 1] is the unit interval. Let C be the middle-third Cantor set. Then 
the corresponding indicator function ipc := xc satisfies the scaling identity (see 
03), ¥>c(f ) = <Pc(x) +^c{x- 2). 

In |45j we use this as the starting point for a multiresolution construction in 
a separable Hilbert space built with Hausdorff measure, and we identify the two 
mother functions which define the associated wavelet ONB. 

Since both constructions, the first one for the Lebesgue measure, and the 
second one for the Hausdorff version (dx) s , arise from scaling and subdivision, it 
seems reasonable to expect multiresolution wavelets also in Hilbert spaces con- 
structed on the scaled Hausdorff measures Tt s which are basic for the kind of 
iterated function systems which give Cantor constructions built on scaling and 
translations by lattices. 
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